Exact dynamics of a reaction-diffusion model with spatially alternating rates 
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We present the exact solution for the full dynamics of a nonequilibrium spin chain and its dual 
reaction-diffusion model, for arbitrary initial conditions. The spin chain is driven out of equilibrium 
by coupling alternating spins to two thermal baths at different temperatures. In the reaction- 
diffusion model, this translates into spatially alternating rates for particle creation and annihilation, 
and even negative "temperatures" have a perfectly natural interpretation. Observables of interest 
include the magnetization, the particle density, and all correlation functions for both models. Two 
generic types of time-dependence are found: if both temperatures are positive, the magnetization, 
density and correlation functions decay exponentially to their steady-state values. In contrast, if 
one of the temperatures is negative, damped oscillations are observed in all quantities. They can 
be traced to a subtle competition of pair creation and annihilation on the two sublattices. We 
comment on the limitations of mean-field theory and propose an experimental realization of our 
model in certain conjugated polymers and linear chain compounds. 

PACS numbers: 02.50.-r, 75.10.-b, 05.50.+q, 05.70.Ln 



I. INTRODUCTION 

Nonequilibrium many-body systems abound in the 
physical and hfe sciences and have recently received much 
attention (see e.g. [3, 12,3 3] and references therein). De- 
spite these efforts, a comprehensive theoretical frame- 
work is still lacking: As yet, there is no equivalent of 
Gibbs ensemble theory for nonequilibrium systems. As a 
consequence, in contrast to equilibrium statistical me- 
chanics, macroscopic observables cannot be computed 
without explicit reference to the imposed dynamics, gen- 
erally described by a master equation, and most progress 
in the field is made by studying paradigmatic models 
0- In this context, exact solutions of simple models are 
scarce, but very precious, since they can serve as test- 
ing grounds for approximate and/or numerical schemes 
and shed light on general properties of whole classes of 
related models. Not surprisingly, nontrivial solutions are 
almost entirely restricted to one dimension (ID; see e.g. 
0, 13 ) , and have focused on completely uniform lattices 
with site-independent rates. Clearly, however, one would 
like to take into account more complex situations, e.g., 
those with spatially varying coupling constants or rates. 
Arguably, one of the simplest generalizations beyond a 
completely uniform system is one with alternating rates. 
In the following, we consider a ID kinetic Ising chain 
(KISC), coupled to two alternating temperatures and en- 
dowed with Glauber-like dynamics. Our analysis of this 
model provides a full description of its dual counterpart, 
namely a reaction-diffusion system (RDS) , characterized 
by spatially alternating annihilation and creation rates. 
Members of these two classes - i.e., kinetic Ising and 
reaction-diffusion models - are prototypical nonequilib- 
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rium systems which have been thoroughly studied on ho- 
mogeneous lattices 0, H S H H • Yet, they still offer 
surprises and novel behaviors, when non-trivial spatial 
rates are investigated. 

Our model was first introduced by Racz and Zia Q 
who recognized that (stationary) two-point correlation 
functions are easily found exactly, even though spins on 
alternating sites are coupled to different temperatures. 
Schmittmann and Schmiiser subsequently realized that 
all stationary correlation functions are exactly calculable 
0. While this information is equivalent to the full sta- 
tionary solution, its representation as exp(— Tioff) is non- 
trivial, involviiig a proliferation of longer-ranged multi- 
spin couplings 10] . Finally, we recently reported the ex- 
act solution for all dynamic correlation functions, start- 
ing from a very simple initial condition, i.e., zero magne- 
tization and vanishing correlations [TH . 

In this article, we complete these earlier studies by 
demonstrating how competing site-dependent rates may 
dramatically affect the dynamics by giving rise to an 
oscillatory approach toward the nonequilibrium steady 
state. We use a generating functional approach to ob- 
tain the complete solution for all correlation functions 
with arbitrary initial conditions. We focus specifically 
on the dynamical magnetization and the spin-spin corre- 
lations and explore their long-time behavior. We will also 
consider the dynamics of domain walls in the spin chain 
which can be mapped onto a reaction-diffusion system. 
Interpreting our results in the language of particle anni- 
hilation and creation, negative "temperatures" acquire a 
natural physical meaning, leading to unexpected oscilla- 
tory dynamics. From a more technical point of view, we 
are able to obtain a complete solution for two nontrivial 
nonequilibrium many-body systems which provides some 
insight into the solvability of two whole classes of related 
models. 

The mapping to a reaction-diffusion system is of inter- 
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est for two reasons. On the theoretical side, the equations 
for densities and correlation functions in the RDS form 
an infinite hierarchy whose solution is not at all appar- 
ent until one recognizes the equivalent spin chain model. 
Also, from an experimental perspective, it is well known 
that diffusion-limited reactions with annihilation and cre- 
ation of pairs of particles are good models for the pho- 
togrowth properties of excited states (solitons/antisoliton 
pairs) in certain coiuugated polymers and linear chain 
compounds 0, 0, Il4j . We propose that spatially al- 
ternating creation/annihilation rates in these systems - 
especially in MX chain compounds - can be generated 
with the help of a laser with spatially modulated power 
output. 

This article is organized as follows: In the next section 
we introduce the kinetic spin chain and its dual reaction- 
diffusion system. Section III presents the complete so- 
lution of the spin chain. Some technical details are rel- 
egated to two Appendices. In Section IV, we map the 
two-temperature spin chain onto a reaction-diffusion sys- 
tem with alternating rates, whose density and correlation 
functions are computed. We analyze the conditions un- 
der which damped oscillations characterize the approach 
to the steady state, and we compare our exact results 
to a simple mean-field description. Section V is devoted 
to a brief discussion of the solvability of related models, 
with Section VI reserved for our conclusions. 



II. THE MODELS 

We consider two closely related nonequilibrium many- 
particle systems on a one-dimensional lattice: (?) a ki- 
netic Ising spin chain (KISC) endowed with a generalized 
Glauber-like dynamics; and (ii) a reaction-diffusion sys- 
tem (RDS), with spatially periodic pair annihilation and 
creation rates. For convenience, we restrict ourselves to a 
periodic lattice (a ring) with an even number of sites and 
study the thermodynamic limit. We expect our exact re- 
sults to be valid for the general cases of odd number of 
sites and/or arbitrary boundary conditions, apart from 
the usual caveats. 

Since the RDS follows from the spin chain via a duality 
relationship, we focus mainly on the detailed description 
of model (i). A spin variable, aj = ±1, denotes the value 
of the spin at site j, with j = 1,2, ...L, and L an even 
integer. Nearest- neighbor spins interact according to the 
usual Ising Hamiltonian: H = —J^jCrj<Jj+i, where 
J > (J < 0) is the (anti-) ferromagnetic exchange 
coupling. Our model is endowed with a nonequilibrium 
generalization of the usual Glauber |^ dynamics: spins 
on even and odd sites experience different temperatures, 

and To, which implies the violation of detailed balance 
To be specific, a configuration {cti, (72, . . . , ctl} 
evolves into a new one by random sequential spin flips: 
A spin (Tj flips to —Uj with rate 

Wj{{a}) = Wj{aj ^ -cTj) 



1 7,- 

= 2 " ^cTj i'^j-i + (1) 

where 72^ = 7e = tanh (2 J/ZcbTg) and 721+1 =70 = 
tanh (2 J/fcfjTo), on even (j — 2i) and odd (j = 2i + 
1) sites. The time-dependent probability distribution 
P{{a},t) obeys the master equation: 

dtP{{<j},t)^ 

= EK(W^')^(W^^)-^.(W)^'(W,^)] (2) 
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where the state {a}^ differs from {a} only by the spin 
flip of aj. Our main goal in this work is to com- 
pute the time-dependent distribution P{{(T},t). To do 
so, we compute all correlation functions (aj-^ . . ■CfjJ)t = 
E{a-} '^h ■ ■ ■ '^in-P({'''}; ^) ^'^d invoke the following rela- 
tionship 0: 

2^P({CT},t) = l + ^a,{a,)t + ^cjjak{cjjak)t + 

i j>k 

+ ^ (Jjak(Jl{(7j(Jk(7l)t + ■■■ (3) 

3>k>l 

This expression illustrates that the knowledge of all 
equal-time correlation functions is equivalent to the com- 
plete knowledge of the distribution function P{{<j},t). 
Recently, this implication was exploited for the steady 
state 9], and for the time-dependent situation yet re- 
stricted to a particularly simple initial condition 

The spin-flip dynamics of this Ising chain can be ex- 
pressed in terms of the creation, annihilation and dif- 
fusion of domain walls, i.e., pairs of spins with oppo- 
site sign. For example, flipping aj in the local conflg- 
uration (Jj-i = ctj = CTj+i = +1 creates two domain 
walls: cTj-i = —CTj and aj = — CTj+i, located on the 
bonds (j — 1, j) and + 1). Similarly, flipping aj in 
the local configuration aj-i = aj = —aj+i = +1 has 
the effect of moving the domain wall on bond (j, j -I- 1) 
by one lattice constant to the left, corresponding to do- 
main wall diffusion. By identifying a domain wall with a 
"particle" , A, our spin-flip dynamics can be recast as a 
reaction-diffusion model, and the two examples translate 
into 00 A A and 0A ^ j40, respectively. The mapping 
from the KISC into its dual RDS is described in detail in 
Table 1. 

Clearly, the presence of alternating temperatures T^, 
To in spin language translates into alternating pair anni- 
hilation and creation rates (1 ± 7e,o)/2 in the RDS. We 
can see easily that letting Tg or To vanish simply prohibits 
pair creation entirely at even or odd sites. Remarkably, 
we can derive an additional, and possibly rather unex- 
pected, benefit from this mapping: Assigning negative 
values for the temperatures Tg and/or To may appear ar- 
tificial in the KISC, but is perfectly natural in the RDS: 
For example, Tg < simply corresponds to a creation 
rate (1 — 7e)/2 > 1/2 which is easily implemented in a 
simulation. In other words, the RDS version is physically 
meaningful, and readily accessible, on a much wider pa- 
rameter space. 
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TABLE L Basic processes underlying the KISC (left) and RDS (middle) dynamics 



Spin flip of site j 


Reactions at bonds next to site j 


Rates 


H » + + -and--H >- + + 


A% — > %A and %A — * A% 


1/2 


H \ > + + + (j oven) 


AA — > 00 (j even) 


(l+7e)/2 


H \ ► + + + (j odd) 


AA — > 00 (j odd) 


(l+7c)/2 


+ + H » H h (j even) 


00 — > AA (j even) 


(l-7e)/2 


+ + + — > + -+ {j odd) 


00 — > (j odd) 


(l-7o)/2 



III. COMPLETE SOLUTION OF THE KINETIC 
SPIN CHAIN 



on the distance between the two sites and their parity, 
^i{k),^i{])&{e,o} [13: 



In this section, we completely solve the dynamics of 
the KISC. It was shown previously that the gener- 
ating function, and hence the full distribution P[{a},t), 
of a broad class of Ising models can be computed from 
two very basic observables, namely: (i) the magnetiza- 
tion, mj{t) = {(Jj)t for arbitrary initial condition, and 
(zi) a particular two-point equal-time correlation func- 
tion, Cj^k{t) = {aj(Jk)t, the resultant from the special 
initial conditions: mj{0) = Cj,fc(0) = (see Appendix A 
for a more detailed discussion of this statement). Here, 
(•)t = J2{a-} ' denotes the usual configurational 

average. In the following, we assemble the necessary in- 
formation about these two observables. 



where 
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V7j7fe (fc - j) 



dr 



l=he+ 7o)/2. 



e-^4_,(ar) (5) 



(6) 



For long times, these settle into their stationary values 
independent of initial conditions: 



CTfCTfe 



Cfe J (oo) 
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V7j-i7fc-i 



(7) 



A. The general t-dependent magnetization. 

From our earlier work pdj , we recall that the magneti- 
zation rnj(t) = {aj)t of the KISC obeys the equation of 
motion, j^mj{t) = ^ [mj_i(i) -|- mj+i{t)] — mj{t) which 
is easily derived from the master equation, Eqn. 10). As 
shown in the general solution of this linear equa- 
tion takes the form mj{t) = J2k ^3,k (t) mk{0), where 
the "propagator" Mj.fc (t) can be written in term of mod- 
ified Bessel functions of first kind In{t) jl^ : 



Af,-fc (t) = e-\ -^Ik-j{(^t), with a = (7e7o)'/' (4) 
V 7fc 

If 7e7o < 0, the propagator becomes Mj,k (t) = 
^(_l)(fc-j)/2|^^./^^|i/2g-tj^_^.(|„|i) where J„(0 is 

a Bessel function of the first kind, with damped oscilla- 
tory asymptotic behavior Tg^. This translates into an 
oscillatory decay of the magnetization . 



B. A special two-point equal-time correlation 
function. 

The second fundamental quantity, i.e., the equal-time 
spin-spin correlation function Ck,j{t), with k > j, is al- 
ready known from For our purposes, it suffices to 
consider an initial condition with zero magnetization and 
zero initial correlations. With the boundary condition 
(cjffc)t = 1 for j = fc, this basic correlation depends only 



where 



a 



1 + 



(8) 



a quantity that reduces to the familiar tanh (J/khT) in 
the equilibrium Ising chain. The approach to these val- 
ues is exponential and monotonic, as e~^'^^~")*t~'^/^, pro- 
vided 7e7o > 0. However, for 7e7o < 0. the approach 
is oscillatory and damped by e~^*t~'^/^ For later 

reference, it is convenient to display the parity depen- 
dence explicitly. Since translation invariance ensures 
^k~i (0 = (^)i W6 need to distinguish three types 



of correlations. The simplest display, which manifestly 
shows the underlying symmetries, is 





/ 7/70 








V 7/7e 



(k-j) 







—e Ik-j (ar). 



(9) 

Note that the last factor is of exactly the same form as 
in the ordinary Ising chain coupled to a single thermal 
bath, the only difference being the geometric mean of the 
two 7's here plays the role of 7 = tanh {2J/ki,T). Before 
turning to the general case, let us remind the reader that 
Eqns. |(SJ) and Q give the time-dependent correlations 
only for a system with no initial magnetization and two- 
spin correlations (e.g., a random distribution). In partic- 
ular, these forms, also used in the next sections, should 
not be confused with the more general cases considered 
in Appendix B. 
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Generating function and general multi-spin 
correlations. 



+ 



In this section, starting from our knowledge of rrij (t) 
and Ckj(t), we c omp ute the generating function of the 
KISC, following [l^- By construction, this generating 
function allows us to find all correlation functions, sub- 
ject to arbitrary initial conditions. A few additional tech- 
nical details are provided in Appendix A. 

The generating function is defined via 'ii{{r]},t) = 
(Yij ~^ Vj'^j))t, where the {r/j} are standard Grass- 
mann variables [l5l fl^ . In the thermodynamic limit, 
L — !■ cx), it simplifies to 



exp I J2 '?Ji^J2 

\j2>jl 



^32,3 



(10) 



If the initial magnetization and all initial correlations 
vanish, the average (...)o on the right hand side of Eqn. 
(|10|l reduces to unity, and one recovers the bilinear form 
for 'i/{{r]},t) which we already reported in Eqn. 
(fTn|l is one of the key results of this paper. 

Given the generating function, all correlation func- 
tions can be obtained by simple differentiation |^ : 

As an illustration, we 



{v}=o 



compute the equal-time spin-spin correlation functions, 
for k > j: 



{crj<7k)t = 



drjkdrjj 



{r,}=0 



= Ck.j (t) + 



(11) 



We emphasize that this is a completely general result, 
valid for any initial conditions, whether homogeneous or 
inhomogeneous, translationally invariant or not. The two 
terms in Hll|l have simple interpretations. While the sec- 
ond term reflects the decay of the initial correlations, the 
first provides the buildup to the final stationary values 
given above iQ. Thus, we see explicitly how the station- 
ary spin-spin correlation function becomes independent 
of the initial values. 

Higher order correlations are can also be evalu- 
ated but are rather complex for general initial condi- 
tions. For uncorrelated, non-magnetized initial condi- 
tions, however, they simplify significantly |llj |. For ex- 
ample, the 4-point function {<Jji<Jj2'^j3crj^)t factorizes 
into two-point functions, according to {o'jiO'j2aj^aj^)t = 
Cj2ji(0cj4j3(^) - Cj3ji(t)c_i^2(t) + Cj4ji(t)cj3j2(i) for 



j4 ^ 33 ^ j2 ^ ji [llj . Similar factorizations 
hold for all correlations. Their steady-state behav- 
ior can be computed directly from the master equa- 
tion or from the stationary limit of the generat- 
ing function, ^{{ri}, oo) = exp ^X]fc>j ^j^fcCfej(oo)y 

Thanks to this simple form, the 2n-point correla- 
tions factorize into a product of 2-point correlations: 



-l"j2n/00 



^31 "32/oc 



^ J2n-l" 32n /OO, 



where j2« > J2n-i > ■ ■ ■ > ]2> ji- 

Finally, following Refs Hill, we can also derive the 
unequal-tmic spin-spin correlation functions Ck,j {f; t) de- 
scribing how a spin on site k at time t is correlated with 
the spin on site j at a later time t + t': 



J 



i 

= J2MjAt')ckAt) Ki'^fi>oA^j,«(i') [Mk,, kit) Ml, At) - Mk,At)Mi,Mt)] (12) 

I 



As an illustration of these general results, in Appendix 
B we specifically compute the spin-spin correlation func- 
tions for general translationally invariant initial condi- 
tions. 



IV. CONSEQUENCES FOR A 
REACTION-DIFFUSION MODEL WITH 
ALTERNATING RATES 

In this section, our exact results will be translated 
into the language of the corresponding reaction-diffusion 
model. We first associate a site j on the dual lattice 



with every bond [j — l,j) of the original chain. Since 
the particles of the RDS are identified with domain walls 
in the spin chain, they obviously reside on the dual lat- 
tice. Each site j can be occupied by at most one parti- 
cle, described by an occupation variable rij which takes 
the value (1) if the site is empty (occupied). Since 
a domain wall involves two neighboring spins, the map- 
ping from spin to particle language is nonlinear, namely, 
Uj — 5 [1 — (Tj^iaj]. As before, we seek the probabil- 
ity, P({n},t), to find configuration {n} at time t, and 
its averages: the local particle density pj{t) = {nj)^ = 
S|n> "j^({'^}j'') s-iid the TO-point correlation functions. 
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("ii ■ • • nj„,)t = '"■h ■ ■ ■ nj„,Pi{n}, t). To simplify 

notation, we continue to denote averages by (•)( for both 
spins and occupation variables, even though they are 
controlled by different statistical weights, P{{a},t) and 
P{{n},t), respectively. In each case, it should be per- 
fectly clear from the context which distribution is rel- 
evant. The dynamics of our model is characterized by 
symmetric diffusion of particles (with rate 1/2) and pair 
annihilation/creation of particles with spatially alternat- 
ing rates (1 ± 7j)/2. In this case, the two particles are 
created on the (dual lattice) sites j and j + 1, by flip- 
ping a spin on the (original lattice) site j. Since jj can 
be positive or negative, subject only to —1 < 7j < 1 
for all j , two very distinct behaviors emerge: (i) when 
both 7e and jo are positive (corresponding to positive 
"temperatures" in the spin model), the annihilation pro- 
cess always occurs with a larger rate than the creation 
process, irrespective of whether j is even or odd; (ii) 
when, e.g., jo is negative and 7e positive, the system dis- 
plays a mild site-dependent frustration: at even sites j 
(i.e., j even and j -I- 1 odd) annihilation is more likely 
than creation, whereas the situation is reversed on odd 
sites (where j odd and j+l even) . As we will see shortly, 
this gives rise to oscillatory dynamics. 



Before diving into the details, some further remarks 
on physical realizations of this model are in order. When 
the rates are uniform (7e = 7o), it is well known that 
such an RDS describes the dynamics of photo-excited 
solitons in conjugated polymers or linear chain com- 
pounds. MX chain compounds, [Pt(en)2] [Pt(en)Cl2]l4, 
where Y stands for CIO4 or BF4 and (en) for enthylene- 
diamine, are of particular experimental interest (12. il3i|. 
In these compounds, photogenerated solitons are so long- 
lived that they can be experimentally studied. Irradia- 
tion with continuous wave (non-pulsed) blue light gener- 
ates soliton-antisoliton pairs which can diffuse apart or 
annihilate. Their static and dynamic properties are in 
quantitative agreement with theoretical models 0, ^3 . 
Since creation, annihilation, and hopping rates can be 
controlled by tuning the laser power, we believe that spa- 
tially alternating rates such as ours will be generated if an 
MX chain compound is exposed to a spatially modulated 
light intensity. 



Returning to our model, our goal in this section is first, 
to derive all correlation functions from our exact solution 
of the KISC. We will also comment on the validity of a 
simple mean-field theory which is widely used for the ho- 
mogeneous (7e — 7o) case JJ, 19J. Further, we show that 
particle hops in the RDS develop a peculiar directional 
preference in the steady state, even though there is no 
explicit bias in the rates, boundary or initial conditions. 
Finally, we illustrate how oscillatory behaviors may result 
from a competition of the underlying processes. 



A. Density of particles in the RDS 

The observable of most immediate interest is the aver- 
age density of particles, Pj{t), in the RDS. Its equation of 
motion can be derived easily from the associated master 
equation, resulting in: 

(2 - 7i - 7i-i) + (7i-iPi-i(*) 



+ l]P]+i{i)) - (4 - I3 - 7j-i)Pi(i) 

- 2[7j(njnj_i)t +7j+i("i"i+i)t] (13) 

It is worthwhile noting that this equation is the first 
member of an infinite hierarchy, connecting lower-order 
correlations to higher-order ones. In general, such hi- 
erarchies cannot be solved directly, without recourse to 
crude approximations. Here, the mapping to the spin 
chain develops its full power, allowing us to compute all 
correlation functions for the RDS. 

The mapping from spins to particles implies that 
Pj{t) = (rij) = ^[1 — ("'i-iO'jOJj so that we can just turn 
to Eqn. to read off the answer. To express it fully in 
RDS language, we also need to translate the initial cor- 
relations, ((Tfccrf)o. For k < I and any t (including t — 0), 
we may write {akcri)t = {ak(Jk+i<7k+io_M-JL:^cri-i'^t)t = 

1 |2C1 whence we 



((l-2nfe+i)(l-2n^2)---(l-2n,-))t | 
obtain, for arbitrary initial condition: 



k<e 



((1 



2"fe+i)(l 



2'^fc+2) •■■(!- 2n^-))o 



X [Mkj-i{t)Me,j{t) - Mkjit)Mi^,^i{t)] (14) 

Since the "propagators" Mij{t) decay exponentially as 
t 00, the steady-state density is independent of initial 
conditions and spatially uniform: 



p{co) = pj{co) 



7 



UJ 



(15) 



\/7e7o 

In Appendix B, we explicitly evaluate Eqn. (|14|l for 
a generic but simple initial condition, characterized by 
a uniform, uncorrelated initial distribution of particles, 
with density p(0) . For simplicity, we discuss only its long- 
time limit here, for p{0) = 1/2. We observe two distinct 
kinds of behaviors: 

(i) When 7e7o > 0, the stationary density of particles is 
approached exponentially fast [except when 7e = 7o = 
±1, see (IB15|) ]. with inverse relaxation-time 2(1 — a) , 
and a subdominant power-law prefactor t~^/^: 



Pit) 



2t 



— e '^Ii{aT) 

T 



p(oo) 



7 

" JO 

^-3/2g-2(l-a)t 



(16) 



2V27ra(l - a) 

This long-time behavior is very similar to that found in 
the usual (je = 7o 7^ il) pair diffusion, annihilation, 
and creation process AA ^ 00 ^ ilSj . 



6 



(ii) For 7o7e < 0, we observe a competition between the 
different processes. For example, when — 1 < 7o < and 
< 7e < 1, the annihilation (creation) reaction dom- 
inates on even (odd) sites. As a result, the stationary 
density is reached exponentially fast with damped oscil- 
lations: 



(17) 



p(oo) 



sin (2|a|< - f ) + \a\ cos {2\a\t - f ) 
4(1 + |a|2)iy^f|^ 



For initial densities other than 1/2, as shown in Ap- 
pendix B, only the amplitude, or the subdominant power- 
law prefactor, of the expressions Ijltil I17|l change. Since 
they depend on all parameters of the model, including the 
initial density, the dynamics is manifestly nonuniversal. 



B. Two-point correlation functions of the RDS 

A deeper understanding of the time-dependent spatial 
structures of our RDS is provided by the m-point cor- 
relation functions, (rij^ . . .nj,„)t of such a model. These 
are related to the correlation functions of the dual spin 



chain, via 



'((1 



(1 



f^j,„-iO'j^))t, and are therefore exactly known. It is in- 
teresting to note that the m-point correlation function 
for the RDS is a superposition of all 2n-point correla- 
tion functions for the spin chain, with n = 1, 2, ...,to. In 
the following, we focus on the most directly observable 
correlation, namely, the two-point function. To avoid 
unnecessary technical complications which add little in- 
sight, we specifically consider a system that is initially 
homogeneously half-filled with A particles, without any 
initial correlations: pj(0) = 1/2 and (7ij(0)ri^,(0)) = 1/4, 
for j ^ k. Such an initial configuration corresponds, in 
the KISC picture, to a system with initially neither mag- 
netization nor correlations. In this case, as we showed in 
[ll[. the generating function takes a rather simple bilin- 
ear form which simplifies the spin-spin correlations. 

With this initial condition, both the spin chain and the 
RDS are translationally invariant, modulo period 2. As a 



result, the two-point correlations C 



it) 



between two sites j and k (with k > j) depend only on 
the distance k — j and the parity /i(fc),/i(j) G {e, o} of 
the two sites. We therefore need to distinguish four dis- 
tinct correlation functions: (t), Cf° U), Cf" .(t), and 

C?° it\ By virtue of our mapping to the KISC, these 

are determined by the 2- and 4-point spin correlations 
as explained in Section III.C [from Eqns. © and Hll)|) ]. 
Exploiting translational invariance, the two-point corre- 
lations for the RDS, for k > j, then follow as: 



(1 



(18) 



Now we are ready to discuss our results. First of all, 
we consider a special case, namely, nearest-neighbor cor- 
relations, li k = j + 1, Eqn. H18|) reduces to 

Again, we should emphasize that the quantities cf°(t), 
c'^{t) and d^it) which appear in this section are the spin 
correlations for a particular initial condition (cf. Eqn. 
©), in contrast to the more general correlations com- 
puted in Appendix B. 

It is interesting to note that, generically, Ci'^{t) ^ 
Cf°{t). Of course, after a little thought this becomes 
less surprising, since {njnj^i)t involves the 4-spin cor- 
relation {aj-i<7jajCTj^i)t = ((Tj_iO'j+i)(. So, if j is odd 
(even), both j — 1 and j are even (odd), leading to a 
contribution of cf^it) vs. C2°{t), respectively. 

For the general case, when k and j are not nearest 
neighbors, this difference between Cf° Jt) and Cf^ Jt) 

o k—j^ ' k~j ' 

does not persist. If k is even and j is odd, we find: 



1 

4 
1 

— c 

4 



(i-cr(t))'-[cr,(i)]^ 



fc-j + l(*)Cfc° 



(20) 



and for k odd and J even, one obtains 



(i-cr(o) 



(21) 



Thanks to the simple relation between even-even and 
odd-odd spin correlations, Eqn. 0, the two right-hand 
sides are now identical. 

A similar line of reasoning shows that C^^ it) = 

k—j 

C-" .{t) for arbitrary separation k — j. Invoking the two- 
spin correlations again, we may write 



(22) 



(1 - cntyf ci%{t)ct_^{t) + cr_^+i(t)cr,_i(i) 



In the following, we discuss the consequences of these 
results. We first consider the steady state. Recall- 
ing our previous analysis of the spin correlations, Eqn. 
(|7|). the stationary limit of the density-density corre- 
lations becomes very simple: Provided k — j > 1, we 



find q . v^; - 4 



(oo) = i(l-cr(oo))^ 



p^{oo)- In 

other words, the two-point correlations of non-nearest- 
neighbor sites factorize into one-point functions, inde- 
pendent of parity. This kind of mean-field-like behavior 
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is typical of free fermion systems 0, 0|. However, the 
nonequilibrium nature of this model still imposes its sig- 
nature. Turning to the nearest-neighbor correlations, we 
find that this simple factorization no longer holds - ex- 
cept in the special case where 7e — Jo- More specifically, 
we find 



and 



V 



(23) 



(24) 



Considering, e.g., < 70 < 7e, we find that Cf°(oo) is 
enhanced over the mean-field result while C°'^(oo) is sup- 
pressed. This can be understood easily: Since 70 < 7e 
implies Tg < To, energetically costly spin flips occur more 
frequently on odd sites j, creating a particle pair on the 
nearest- neighbor dual sites (j-l- 1, j). Clearly, these sites 
form an (e, o) pair. Moreover, the rate for pair annihi- 
lation is lower on (e, o) sites. Hence, particle pairs are 
more likely to reside on (e, o) than on (o, e) sites. This 
also implies that (e, o) sites act as net particle sources, 
while (o, e) sites function as sinks |10| . Not surprisingly, 
therefore, we find Cf°{oo) > C5"^(oo). By virtue of this 
reasoning, it is also immediately apparent that this dif- 
ference can only persist for nearest-neighbor correlations. 
The same argument holds for 70 < < 7e. 

A direct consequence of C°'^(i) 7^ Cf°{t) is the presence 
of a peculiar directional preference in the RDS. If we 
consider a particle on site j, we can ask for the average 
rate, TZj(t), with which it will jump to the left (i.e., to 
site j — I) vs to the right, defined as TZj{t) = ^{nj{l — 
nj+i) — nj(l — nj_i))t. Here, the first (second) term is 
the average rate for a particle on site j to jump to site 
j+1 {j— 1). In our case, one might expect this difference 
to vanish since neither bulk rates nor boundaries impose 
a directional bias. Moreover, to avoid a potential bias 
at t = 0, we start from a translationally invariant initial 
condition with p{0) = 1/2 . Yet, since 7^J(^) cx Cl^it) - 
Cl°{t), it is manifestly nonzero. Explicitly, we find: 



7^J(t) 



2a _ 1 

7e 



cf{t) ,ieven 
cf{t) ,iodd. 



(25) 



this expression reduces to a density difference which van- 
ishes for all times provided the initial condition is ho- 
mogeneous. For inhomogeneous initial condition, J7}(t) 
exhibits nonzero transients for finite times but decays as 
t 00. 

Let us conclude this section with a few brief remarks 
about the validity of the mean-field approximation for 
this system. We already noted that it does not predict 
the nearest-neighbor correlations correctly, except in the 
special case je = 7o- We now show that it also generically 
misses the stationary density. 

We begin by recalling Eqn. H13|l . Seeking a trans- 
lationally invariant (modulo 2) solution with P2j{t) = 
Pe{i), P2j+i{t) = Po{t) for all j, the mean-field approx- 
imation corresponds to truncating two-point functions: 
(nj(t)nj±i(i)) ~ pe{t)po{t). Starting from a uniform ini- 
tial density p(0) , we find pdt) = Po(0 = PMF{t), with 



PMF{t) 



Pp [p{0) - Pm] - P,n [p{0) - Pp] e-W^-h^+f^y 

p(0) - prn + [P(0) - Pp] e-V4-(7.+7.)^ 



where 



Pp,': 



2 ^ \/4-(7e+7o)2 



le+lo 



le+lo 



which even persists in the steady state: 



The stationary limit is clearly pj^ip (oo) = pp which dif- 
fers from our exact result, Eqn. H15|l . except if 7e = 7o. 
In other words, the remarkable accuracy 0, |3 of 
the mean-field approximation for the stationary state of 
the uniform system (72 = 70) appears to be an "acci- 
dent" due to the fact that when rates are uniform the 
steady state is a product measure. We also note that 
the exact relaxation time to the steady state, T^xact = 
[2 — ^7e7o]~"'", does not coincide with the mean-field pre- 
diction, tmf = [\/4 — (7e -I- 7o)^]~"'"- For such dynamic 
quantities, the exact and the approximate results differ 
for any choice of 7e and 70. In particular, the mean- 
field theory always predicts an exponential decay to the 
steady state, completely missing the possibility of oscil- 
latory behavior. 



SOLVABILITY AND RELATIONSHIP WITH 
FREE FERMION SYSTEMS 



, j even 



Specifically, for 70 < 7e, particles on an even (odd) site 
jump preferentially to the right (left). Of course, this 
directional preference vanishes as soon as 7e = 7o ■ More- 
over, even when it is nonzero, it does not generate a 
mass current. Counting the net flow of particles between 
sites j and j -I- 1, the natural definition of such a cur- 
rent is Jj{t) = ^((1 — nj)nj^i — nj{\ — nj+i))^. Clearly, 



The crucial ingredient for the solvability of the KISC is 
the quadratic spin dependence of its Glauber-like kinet- 
ics. Thanks to this simple form, the hierarchy of equa- 
tions for the correlation functions is closed: to solve the 
equations for the A^-spin correlation functions, one needs 
to know only m-point correlations with m < n. 

In RDS language, the dynamics of the particles can be 
rewritten as a free fermion model, by defining a suit- 
able quadratic (but non-Hermitian) "stochastic Hamil- 
tonian" . Following standard methods 0, 0, ITsL l2l| , we 
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can rewrite the master equation for the RDS as a formal 
imaginary-time Schrodinger equation: (d/dt)\P{t)) — 
—H\P{t)). The Hamiltonian H is constructed by associ- 



ating the usual Pauli matrices a- {erf) with the creation 
(annihilation) of a particle at site j : 



'2^ = E [4^7+1 + 4+1 + (1 + 7e)44+i + (1 - Ie)a7a7^, 7e(^74 + - (1 



J ei;en 



E + 4+1 + (1 + '7°)44+i + (1 - '7°)'^i ^i+i - 4 + - (1 - 7o) 



(27) 



The key to the solvability of this Schrodinger equa- 
tion lies in the bilinear dependence of the Hamiltonian 
on the Pauli matrices. This is due to the fact that the 
spin-flip rate s itTll implicitly fulfill the free fermion con- 
straint 13, 0, llSf - In RDS language, this condition re- 
quires that the sum of the particle diffusion rates equal 
the sum of the (local) annihilation and creation rates, 
i.e., 1/2-^1/2= (l-h7j)/2-h(l-7j)/2 withj € {e, o} in 
our case. If this relation is violated, H includes quartic 

terms, of the form '^J '^t '^j+T-'^t+'i-' ^^'^ ^^'^ associated 
RDS can no longer be solvecf exactly. It can, of course, 
be simulated, and for those cases investigated so far, it 
appears that the quartic terms are irrelevant for the long- 
time dynamics ^ .IQl . It is also worth noting that the 
free fermion constraint is not particularly artificial: the 
simplest models for photogenerated solitons in MX chain 
compounds satisfy it quite naturally 12]. 

Here we decided to invoke generating function tech- 
niques instead of diagonalizing (|27(l . In our view this is 
the most convenient and systematic approach to solve si- 
multaneously both the KISC and RDS, for two reasons. 
First, the free- fermion approach requires various techni- 
cal steps {e.g. introduction of so-called pseudo-fermion 
operators and a Bogoliubov-like transformation) which 
make the general treatment rather involved, especially 
for arbitrary initial conditions 0,0] . Further, the diag- 
onalization of (|27|) yields only correlation functions with 
an even number of spins (see Section III and IV); the 
calculation of correlations involving an odd number of 
spins requires a dual transformation of H27|l into a new 
stochastic Hamiltonian which must also be diagonalized 

M- 

Let us also mention that damped oscillatory decay has 
been observed before in certain reaction-diffusion mod- 
els However, those models, and hence the physical 
mechanisms leading to the oscillations, are completely 
different from ours. As an example, a diffusion-limited 
fusion model is defined by three processes: (i) bi- 
ased diffusion: with rate D(l + t;), ^ 
with rate D{1 — -q) (with < 77 < D)\ {ii) biased fu- 
sion: AA %A with rate D{l + 2ri), AA A% with 
rate D{1 — 277); and [iii) homogeneous pair production: 
00 AA , with rate D. With this special choice of 
rates, the equation of motion for the density closes and 



becomes solvable. In order to observe oscillatory decay 
of the particle density, the initial condition must be in- 
homogeneous. For a homogeneous initial condition, the 
density decays exponentially. In contrast, the equation 
in our reaction-diffusion model does not close, and the 
oscillatory behavior is generic: it occurs for any initial 
condition, inhomogeneous or not. 

VI. CONCLUSIONS. 

To summarize, we have presented a full exact solution 
for the dynamics of a non-equilibrium Ising spin chain, 
with arbitrary initial condition. The model is character- 
ized by a generalization of Glauber dynamics: spins on 
even/odd sites are coupled to alternating temperatures, 
Tf, and Tq. We obtain all correlation functions from a 
generating functional. As an illustration, we have dis- 
cussed the equal-time and the two-time spin-spin corre- 
lation functions. 

Identifying domain walls in the spin system with parti- 
cles on the dual lattice, the model can also be interpreted 
as a reaction-diffusion system. Particles are created and 
annihilated in pairs; the rates for these processes alter- 
nate from even to odd sites. This mapping opens up an 
interesting extension of parameter space: while negative 
temperatures are unphysical for the spin chain, the corre- 
sponding rates are perfectly natural in the context of the 
RDS. By expressing particle-particle correlations as su- 
perpositions of spin-spin correlation functions, the RDS 
becomes exactly soluble. This is not entirely trivial since 
the BBGKY 22] hierarchy for the RDS is not closed: 
its solution is far from obvious unless one recognizes the 
connection to the spin chain. 

For < 7e7o, energetically favorable spin flips always 
dominate over unfavorable ones, irrespective of whether 
they occur on even or odd sites. In RDS language, pair 
annihilation is always more probable than pair creation. 
As a consequence, we find that all quantities decay ex- 
ponentially to their steady-state values. In contrast, for 
7e7o < 0, we observe (damped) oscillatory behavior. Its 
origin can be traced to a competition of pair creation and 
annihilation on even vs odd sites on the original lattice: 
If, say, 7o < 0, then pair creation dominates over annihi- 



9 



lation on odd sites while the relation is reversed on the 
even sites. Hence, a given initial particle density may 
first decrease, due to annihilation processes and then re- 
cover, as the available empty sites are (partially) filled 
again by the strong creation process, and so on, until the 
stationary density is reached. 

Remarkably, even in the absence of any bias in 
the rates, boundary or initial conditions, particles still 
"know" the difference between right and left: For, e.g., 
7o < 7ej particles on even (odd) sites jump preferentially 
to the right (left). Even though this directional prefer- 
ence does not lead to a systematic particle current, it is 
still somewhat surprising. However, once we recall that 
particles are most often created (annihilated) on pairs of 
neighboring sites, with the odd site on the left (right), 
we recognize that the directional preference is simply a 
response to this density gradient. 



Since exact solutions, especially of a full nonequilib- 
rium dynamics are rare, we hope that our model can 
serve as a testing ground for various generalizations or 
approximations. The features reported here - exponen- 
tial decays, damped oscillations, and directional prefer- 
ence - should be generic for a whole class of genuine 
out-of-equilibrium models. Moreover, they should be ex- 
perimentally observable in MX chain compounds exposed 
to spatially modulated laser light. 
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APPENDIX A: THE DERIVATION OF THE GENERATING FUNCTION 



In this appendix we provide some details for the derivation of the generating function (|10() . which is one of the 
key results of this work. We follow Aliev's work and notation jl^J. Aliev established that the generating functions of 
a very general class of disordered Glauber-Ising spin chains, including our case, can formally be expressed in terms 
of two functions MA,(i) and two additional quantities W^i^{t), which depend in a very involved fashion on M^-(t). 
Below, we will see that these quantities are closely related to physical observables, namely, the magnetization and the 
two-point correlations. Here, we follow Aliev by noting that the Laplace transform of M^^, is the inverse of an L x i 

band matrix (s + 1)1 — . For our case, the entries of this matrix can be taken from the rates and read explicitly 
{L is even): 



+ 1)1- v± 

\2j 


-l.k 




To 

1- l)<52j-l,fe - Y(^2j-l,fc-l + S2j- 


_i,fc+i), (1< j < L/2) 


(s + l)l-l[/± 


2j.k 




1- ^)S2j,k — -y (^2j,fe-l + S2j,k+l ) 


{l<j< L/2) 


(s + l)l-ic/± 


l,k 




y(^2,fcT'5fcx) 




(s + l)l-i{/± 


L,k 




hl)5L,fc-y(^L-i,fcT'5i,fc) 





(Al) 



Given HAljl . it is easy to evaluate the inverse of [(s + 1)1 



^U^kk: 



M, 



k.J 




i{k-j) 



a cos 



(A2) 



where 



7r(2Tt-l) 
L 



and = -T^, with n = 1,2, . . . , L. In the thermodynamic limit L — > oo, the two quantities 



M^i. and Mj j, coincide whence we simply have 



k,j 




Mk., ^ J— TT- 



^i(k-j)il> 



2tt s -|- 1 — a cos ( 



(A3) 



Taking the inverse Laplace transform of (|A3p . we recover Eqn. Q for the propagator. 

Since the W^i^{t) can be expressed in terms of the M^i^{t), we may immediately conclude that W^,^{t) Wj,kit) 
as L — !■ oo. According to Aliev, Wj^kit) is simply the two-point correlation function, {ajcrk^j)t, for a particular initial 
condition, namely, TOj(O) = and {ajcrk^j)o = 0. For our case, these correlations were given in Eqn.jSJ). 
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For readers familiar with Aliev's work [Tjl , these remarks fill in the gaps between Aliev's formal and general analysis 
and the special case we are interested in here. It follows that the generating function of our KISC admits the compact 
and explicit representation of Eqn. (jl()|l which encodes the complete dynamics of the system. 

APPENDIX B: TRANSLATIONALLY-INVARIANT INITIAL CONDITIONS: THE SPIN-SPIN 
CORRELATION FUNCTIONS AND THE PARTICLE DENSITY 

In Section III.B, we have derived an exact expression, Eqn. Hll|) . for the spin-spin correlation functions of our 
KISC, valid for arbitrary initial conditions. Here, we impose a natural restriction, namely, translational invariance, 
on the initial conditions. Thanks to the symmetry, Eqn. Hll|) simplifies considerably, as we will show now. 

As we already pointed out in pj| . for translationally invariant initial conditions, we only need to consider the 
correlations between spins at two even sites, two odd sites, and one even, one odd site. We denote these by C2^(t) = 
{(^2iO-2{£+n))t, cg°(i) = {cr2i-i(T2i-i+2n)t, C2°_i(t) = c^n-ii^) = {o'2icr2i+27i~i) t = {cr2i+io-2i+2n) t- Of course, there is 
no need to study n < cases. For the special case of zero initial magnetization and correlations, these correlations 
are already known pj] | and are given by Eqn. 0. Here, we seek their form in a more general case, starting from a 
homogeneous initial condition. 

Let us recall from |Tl| that in the translationally-invariant case the quantities a2n (t) = 
k {leC2n (t) + loC2n {t)] , a2n-i (t) = Q;C2°_i(t), obcy the foUowing simple equation: -^aj = —2a,- + a[aj_i + 
Uj^i] , J > with the initial condition oq (t) = 7 and 7 = (7e + 70) /2. The equations of motion of the KISC pJ| 

also give the following relationships among the correlators: (i) = -^c^nit) + (c2n{0) — ^C2^(0)^ e^^*. 

The explicit expressions for the correlators follow from (|ll(l . or by the methods of images directly from a„(t) = 
On (00) + e~^* X]m>o[°'™(0) ~ Om(oo)] {/„_m(2Q;t) — /„+rn(2Q;t)} where a„(oo) = 7 w'^ and n > 0. From the definitions 
of a„, we immediately infer: 



lo \lo 

— \ X! ["2m (0) - a2m(oo)] { /2(„_™) (2ai) - /2(n+m)(2a0} 

^° ^° rn>0 

ae 



+ E [c2™-i(0) - c§°„-i(«3)] {/2(„-™)+i(2ai) - /2(„+™)-i(2a<)} 

'° m>0 

(7ec^°(0)-7ocr„(0)). (Bl) 



270 

Following the same steps for c'^^-i — '^2n~iJ obtain, for {n > 0) 



cTn-l{t) - 4n-l{^) + e^'' Y.^Cr^-^{0) - cZ^-li^)] {/2(„-m)(2ai)-/2(„+„_i)(2ai)} 

m>0 

„-2t 

E [a2m(0) - a2m(oo)] {/2(n-m)-i (2ai) - /2(n+m)- 1 (2a<) } . (B2) 



a 

m>0 



The expressions (jBl|) - (IB2p illustrate that the time-dependence of the spin-spin correlation function depends non- 
trivially on the initial condition, and we may therefore anticipate non-universal behavior. Of course, when 7e — jo, 
the expressions (|Bip . (jB2p coincide with those obtained by Glauber 

An interesting situation occurs when, say, 70 is negative while < 7e < 1, so that a — i\a\. Then, we have 
l2n{2i\a\t) = (-l)"J2„(2|a|0 and /2„±i(2i|a|t) = ±i (-1)" J2„±i(2|a|t), where J„(x) = ^ cos {x sin q - nq) is the 
Bessel function of first kind jl^]. Further, when a = i|a|, the expressions (|B1|I . (|B2|I become: 



7o \lo 



e 



,-2t 



'^2f^ - ^ E [«2™(0) - «2™(00)] (-1)"+™ { J2(„_„)(2|a|i) - J2(„+„.)(2|«|t)} 
l"^"! m>0 
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+ E['^2™-i(0)-c^™-i(oo)](-l)"+" {J2(„-™)+i(2|a|0 + J2(„+™)-i(2|a|0} 

1'°! m>0 

- ^ (7ec^°(0)-7oc^:i(0)). (B3) 
cZ^iit) = c-_i((X3) + ^ [c- _i(0) - (-1)("+'") { J2(„_„)(2|a|t) + J2(„+™_i)(2|a|t)} 

m>0 

- ^ E [«2™(0) - a2™(oo)] (-1)("+") { J2(n+m)-i(2|«|i) - J2(„-™)-i (2|a|t) } . (B4) 

' ' m>0 

In the long-time limit, these expressions exhibit a damped oscillatory approach to the stationary state. 

We now turn to the equivalent RDS. As we have seen in Section IV. A , the density of particles is related to the 
nearest-neighbor spin correlations, cf, according to p(t) = i(l — c\°{t)). Here, our goal is to determine the long-time 
behavior of this density for a homogeneous (but otherwise arbitrary) initial concentration of particles p(0). In this 
respect, the expressions (|B2(I and (IB4(I are not very practical as they involve infinite sums of Bessel functions. At this 
point, for further convenience, it is useful to introduce four auxiliary functions defined as follows (with < A < 1): 



F,{\t) = ^ A2'"-ie-2*|/2(„_,)(2at)-/2rn(2ai)} =2A(1 + A2) / ^- 



m>0 



+ A4-2A2 cos2g 



m>0 



(1 — Q COS q) 



sin 2q sin q 



1 + A4 - 2A2 cos2g 



cos^ q 



(B5) 
(B6) 

(B7) 

.l + A4 + 2A2cos2(7j ^^^^ 

To establish these expressions, we have invoked the integral representation of the Bessel functions and the 
properties of geometric series. With these functions and the help of Eqn. (|14|l . the density of particles in the RDS 
model can now be recast in compact form. Two cases emerge naturally: 



G2{X,t) 



- ^ A^™-! (-l)-e-2* { J2(i_„.)(2|a|i) + J2™(2|a|i)} 

■m>0 

2A(1 + A2) r^e-^' cos(2|a|ising) 
Jo 7^ 



1 + A4-I-2A2 cos2g 
E (-l)"e-2* {J2™+i(2|a|i) - ./i_2m(2|«|i)} 



m>0 



/ sin (2|a|t sing) 

Jo 



sin 2q cos q 



• When 7e7o > : 

p{t) - p{oo) 

• When 7e7o < : 

p{t) - p{oo) = 



[Fi{Lu,t)+F2iLu,t)] - i [Fi(1 - 2p(0),i) + ^i^2(l - 2p{0),t) 
2a 2 1 a 



-^[G2(l-2p(0),t)-fGi(u;,0]-i 
2\a\ 2 



-G2(w,i)+Gi(l-2p(0),t) 



The case 7e7o = is special and gives rise to a purely exponential time-dependence: 



Pit) 



7 



7-2 



p(0) e 



-2t 



(B9) 



(BIO) 



(Bll) 



We now proceed with the analysis of the long-time behavior of these expressions. Again, we first consider the case 
where 7e7o > and then 7e7o < 0. 

• When 7e7o > 0, the main contribution to the long-time behavior arises from the small q contribution in the 
expression of the functions Fi and F2. Therefore, one may expand the integrand of Fi and F2 in Eqn. I|B9|) . It 
is also essential to pay due attention to the initial condition: 
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1. When < 7e7o < 1 and < p(0) < 1 (also p(0) ^ p{oo)), we obtain 

7^ (l + i){l-2p(0)} + 2p(0)2 



p(t) - p(oo) 



a(l -cj2) 

2. When < jelo < 1 and p(0) = 0, we find 



p(0)(l-p(0)) 



g-2(l-a)t 



pit) - p(oo) ~ - 



1 



7 
a 



3. When < 7e7o < 1 with p(0) = 1, we have 



p{t) - p{co) 



1 



1 

a 



a(l -cj2) 



1 

at 



1 

at 



,-2(l-Q)t 

Ay/nat 



,-2(l-Q)t 

4-\/7rai 



(B12) 



(B13) 



(B14) 



These results show that for 7e7o > 0, the density generically approaches its stationary value as 
cx t~'^/2 g-2(i-Q)t^ with some nontrivial amplitude. Only if the lattice is initially completely empty/occupied 
by particles, the long-time behavior is modified to oc t"^/^ g-2(i-a)t (pj-Qvided ^ ^ ±1)- 
4. The case where 7e_= 7o = 7 = ±1 is critical and we can check from ljB9|) that one recovers the previously 



known results 



pit) - p(oo) 



7 



2V7rt 



(B15) 



In this case, it is well known 0,^3 that the density of particles approaches the steady state algebraically 
slowly (cx t~^/^). For 7=1 (only pair annihilation), the stationary value is p(oo) = 0; in contrast, we find 
p(oo) = 1 for 7 = — 1 (only pair creation). We emphasize that for such a critical dynamics neither the 
dynamical exponent nor the amplitude of (|B15|) depend on the initial condition. 



• When 7e7o < 0, it is difficult to directly analyze the long-time behavior of the oscillating function Gi and G2', 
instead, we seek upper and lower bounds for A 7^ ±1. We observe that the denominator of the integrand in the 
expressions for Gi and G2 can be bounded as follows: (1 - A^)^ < (H- A"* -f 2X^ cos2g) < (1 -I- A^)^ Therefore, 
we obtain for the auxiliary functions Gi and G2: 



A 



-2t 



Ji{2\a\t) 



1 + A2 2\a\t 
' J2{2\a\t) 



< Gi(A,i) < 



A(l + A2) e~2* Ji(2|a|i) 



2\a\t 



1 + A2 



2|a|i 



< G2(A,t) < 2 



A 



1-A2 



2 „-2t 



e-' KM2\a\t) 
2\a\t 



(B16) 

(B17) 



If A = ±1, one has the exact expressions: Gi(l,i) = — Gi(— l,t) 
Ji(2|a|f). At long times and for finite n, e~2* j„(2|Q;|i) 



e-2* Jo(2|a|0; G2(l,t) - G2(-l,i) = 
cos (2|Q;|i — f (2n + 1)) and therefore 



the upper and lower bounds in Eqns. (|B16|I and (|B17(I display the same time-dependence. With the help of 
Ean. (|B10|l . we thus deduce: 



(B18) 



1. When < p(0) < 1 (and obviously p{0) ^ p(oo)), 

p{t) ~ p{oo) ~ t-3/2e-2t (acos (2\a\t + ^) + Bcos (2\a\t ~ 



where A and B are some amplitudes depending nontrivially on all the parameters of the system and on 
the initial density. 

2. When p(0) — 1, we obtain an explicit expression for the long-time behavior of the density: 



p{t) - p{(X)) 



-M2\a\t) 



cos 



(2H^-l)- 



2 2y/^t 

3. When p(0) = 0, we also have an explicit expression for the long-time behavior of the density: 



p(i)-p(oo)~ ^Jo(2|a|t) 



2y^TT\a\t 



I „i I 37r 
cos ( 2\a\t + — 



(B19) 



(B20) 
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These results show that, for 7e7o < 0, the density displays oscillations which are damped by a factor i~''e~^*, 
where j3 = 3/2 for generic initial densities p(0), with two exceptions: we have /? = 1/2 if the system is initially 
completely empty or occupied. 

Of course, following the same approach one would bo able to compute every n— point correlation functions, for 
both, the KISC and RDS. While perfectly straightforward, these computations become rather tedious. 
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